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Abstract 



We consider the problem of critical gravitational collapse of a scalar field 
in 2 + 1 dimensions with spherical (circular) symmetry. After surveying all 
the analytic, continuously self-similar solutions and considering their global 
structure, we examine their perturbations with the intent of understanding 
which are the critical solutions with a single unstable mode. The critical 
solution which we find is the one which agrees most closely with that found 
in numerical evolutions. However, the critical exponent which we find does 
not seem to agree with the numerical result. 
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I. INTRODUCTION 



Gravitational collapse is an inherently dynamical and multi-dimensional problem. Its 
complexity is indirectly revealed by both the relatively few analytic solutions available for 
study and the challenges associated with the construction of a good numerical code to sim- 
ulate the collapse process. Nonetheless, work over the last decade has revealed surprisingly 
rich behavior in even nominally simple (e.g. spherically symmetric) systems Q. These 
efforts have led to important advances in our understanding of the process of black hole 
formation and the presence of "critical" behavior in these dynamical, gravitating systems. 
Indeed, this behavior has been shown to be present in a great variety of systems 0. 

Some more recent work, from which we take our cue in this paper, has been the con- 
sideration of the gravitational collapse of a minimally coupled scalar field in the presence 
of a cosmological constant but in a lower dimensional spacetime, namely 2 + 1. There are 
several motivations for studying such a model beyond the intrinsic interest in examining 
critical behavior in another system. Among these is the recent flurry of work on anti-de 
Sitter (AdS) spacetimes stemming from the AdS/CFT conjecture. This conjecture posits a 
correspondance between the gravitational physics in an AdS spacetime and the physics of a 
conformal field theory on the boundary of AdS. Hence, understanding AdS spacetimes can 
potentially yield insight into Super- Yang-Mills theory (and vice versa). Another motivation 
for studying 2 + 1 scalar field collapse is partly the relative simplification that results in 
going from 3 + 1 to 2 + 1 dimensional gravity. By itself, this wouldn't necessarily be that 
compelling, but there are, of course, some intriguing solutions in 2 + 1 such as the BTZ black 
hole that closely parallel the black hole solutions of 3 + 1 gravity. Earlier work has consid- 
ered the question of gravitational collapse to a BTZ black hole, but using either null fluid or 
dust as the collapsing matter Thus, considering this particular model of gravitational 
collapse could potentially pull together several distinct ideas. 

The study of this particular system has included numerical simulations of the spherically 
(or circularly) symmetric model as well as analytic investigations. The numerical results 
HU have demonstrated the existence of critical behavior at the threshold of (BTZ) black 
hole formation, found a continuously self-similar (CSS) critical solution and calculated the 
value of the mass-scaling exponent.^ On the analytic side, a family of exact solutions to 
the dynamical equations was found under the explicit assumption of self-similarity (and 
zero cosmological constant) 0. In addition, comparing these solutions to the numerical 
results provided strong evidence that one of the solutions in this family corresponded to 
the critical solution found in the simulations. This exact solution is interesting, of course, 
not least because it would appear to be the first analytic solution which has been shown 
to correspond with a critical solution found from numerical simulations. There is a CSS 
solution for the gravitational collapse of a spherically symmetric scalar field in 4D but it is 
known not to be a critical solution ||. 

Our efforts here are focused on the exact, CSS solutions in this model and on identifying 
those that are relevant to critical collapse. The relative tractability of 2 + 1 gravity allows 
us greater analytic control of the dynamics in this system which will hopefully yield to 



In H the exponent was found to be 7 ~ 1.2 while in || it was found to be 7 ~ 0.81. 
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some added insights concerning the collapse. Indeed, as we will see, even the equations 
associated with the perturbations of the CSS solutions are integrable. As a result, we are 
able to attempt a fairly exhaustive study of the CSS solutions and their perturbations. We 
begin in the next section with the general equations of motion for our model. Section III 
examines each of the three classes of CSS solutions and their global structure. In section IV 
we perturb the exact solutions and examine the stability properties of each of these three 
classes of solutions. We also calculate the mode structure of these solutions and identify the 
particular exact solution which has a single unstable mode. Not surprisingly, this turns out 
to correspond to the same critical solution which Garfmkle found most closely resembled the 
critical solution that appeared in the numerical simulations J7|. Sections V and VI consider 
briefly the implications for scalar field collapse when there is a potential for the scalar field 
and offer some conclusions. An Appendix makes some comments on some more general 
solutions and properties of our system in 2+1. 

In concluding this section, we mention that as our work was nearing completion, Garfmkle 
and Gundlach [[| presented some recent perturbation results in this same system. They, too, 
determine analytically the mode structure of these exact solutions with an eye to calculating 
the mass-scaling exponent and confirming that the critical solution has a single unstable 
mode. As such, there is overlap between their work and what we present here. 



II. THE EINSTEIN-SCALAR FIELD EQUATIONS 

The action for a scalar field minimally coupled to 2 + 1 gravity is given by0 



V a 0V> + 2V{<f>) + 



A 



AttG 



(1) 



where R is the Ricci curvature scalar, V(4>) is the potential of the scalar field 0, and G and 
A are respectively the gravitational and cosmological constants. The covariant derivative is 
denoted by V a . The gravitational field is then governed by the Einstein field equations, 



Rab — ~9abR + Ag a 
with the energy-momentum tensor given by 



8nGT ab , 



(2) 



T ab = V a 0V fe - -g ab [V c 0V c + 2V ((/>)] . 
The evolution of the scalar field is determined by 



(3) 



(4) 



2 We will choose units such that the speed of light is one. 
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FIG. 1. The spacetime in the (u, «)-plane. Each of the four regions is defined as follows: 
I = {x a : u < 0, v > 0}, II = {x a : u < 0, v < 0}, I' = {x a : u > 0, v < 0} and 
II' = {x a : u > 0, v > 0}. The null geodesies, l a , lying along the hypersurfaces u = const, 
are outgoing, while the null geodesies, n a , lie along the hypersurfaces v = const, and are ingoing. 

The general metric for a (2 + l)-dimensional spacetime with spherical symmetry^ can be 
cast in the form, 



where (u,v) is a pair of null coordinates varying in the range (—00,00), and 9 is the usual 
angular coordinate with the hypersurfaces = 0, 2tt being identified. It should be noted 
that the form of the metric is unchanged under the coordinate transformations, 



In the following we will use this gauge freedom to fix some of the integration constants which 
arise. In addition, the roles of u and v can be interchanged. To fix this particular freedom, we 
choose coordinates such that along the lines of constant u the radial coordinate r increases 
towards the future, while along the lines of constant v the coordinate r decreases towards 
the future. This, of course, just defines u as outgoing and v as ingoing null coordinates [cf. 




(5) 





(6) 



Fig.l]. 



Defining null vectors, l a and n a , along each of the two rays by 



3 Sometimes the spacetimes described by metric (||) are also referred to as having axial or circular 
symmetry. 
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l a = 6 u a , n a = 5l (7) 

we find that 

l b V b l a = = n b V b n a , (8) 

i.e. these null rays are afnnely parameterized null geodesies. The expansion for each is given 
by 

B l = -V A / A = e - 2CT ^, 
r 

9 n = -V x n x = e' 2 ^, (9) 
r 

where, as usual, the comma notation denotes partial differentiation. 

The corresponding Einstein-scalar field equations (0) for the metric (f|) take the form, 

r )UU - 2a n r u = -8nGr<fr 2 u , (10) 
r >vv - 2a >v r >v = -8irGr(f> 2 v , (11) 



where 



+ 2ra uv = -87rGV {<j> tU< j) <v - 2e 2a v) , (12) 
r> = 87iGre 2a V, (13) 



V{cj>) = V{cj>) + (14) 
The equation of motion for the scalar field (^) is now 

20,™ + - {r iU <t>, v + r^,„) = -e 2a ^p-. (15) 
r o(p 

Because of the presence of the effective potential V(4>), the above equations do not, in 
general, allow the existence of exactly self-similar solutions. However, because our main 
interest is in phenomena occuring at the threshold of black hole formation, we may still seek 



for solutions that are asymptotically self-similar (H|. Indeed, what we are really looking 
for in the study of (Type II) critical collapse is a solution that approaches self-similarity as 
the collapse is about to form a spacetime singularity, that is, in the limit of small spacetime 
scale. To see the above clearly, let us first introduce the dimensionless variables, x and r, 
via the relations 

where uq is a dimensionful constant, and the above relations are assumed to be valid only 
in the region v > 0, u < 0. We will refer to this region as Region I (cf. Fig 1.). Later, we 
will consider extensions of these coordinates to other regions of the (u, v )-plane. In terms of 
x and r, the Einstein field equations (JTSHrSI) take the form 
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"5, xx ~\~ 2S tTX -|- S >TT S jX S,t 2{jJ x -|- CT l r)('S J x S,t 

SirGs (<p jX + <p iT - cf (17) 
s jXX - (1 + 2a jX )s )X = -8ttGs<p% (18) 

+ s tTX - s >x + 2s (a >xx + a )TX ) = SttGs <p >x {y tX + tp )T - c) 

-2u 2 e- 2T+2a+x V((f))] (19) 

s, xx + s, TX - s tX = SttGsuI e - 2T+2 ° +x V(<j>) (20) 
2s ((p jXX + tp jTX ) + tp jX (2s jX + s >T - s) = -tp jT s, x + cs :X 

su 2 e- 2 ^^, (21) 



where 



r(u, v) = —us(u, v), 

(j){u, v) = cln |— u\ + ip(u, v), (22) 

with c being an arbitrary constant. 

In every example of Type II critical collapse studied so far, the critical solutions have 
been found to have either discrete self-similarity (DSS) or continuous self-similarity (CSS) 
0. Discrete self-similarity of the spacetime corresponds to periodicity of solutions in the 
coordinate r, 

A(t,x) = A(A + t,x), (23) 

where the vector A is defined as A = {a, s, <p}, and the constant A denotes the period of 
the solutions, while continuous self-similarity corresponds to the case in which 

A(r,x) = A(x), (24) 

with the homothetic Killing vector £ a being given by 

e = < + <• (25) 

In our current problem, we can see from Eqs. (P"7|-|2"TD that the field equations contain an 
explicit factor of e T and consequently neither of the two kinds, DSS and CSS, of self-similar 
solutions is allowed. The physical reason, of course, is that the presence of the effective 
potential introduces a scale into the problem and, as a result, the field equations are not 
scale invariant. However, as mentioned previously, critical collapse is only relevant in the 
region where r ^> 1 (0 < — u <1). Therefore, provided that 

e- 2T V(cj>) ^ 0, (26) 

as r — ► oo, solutions that asymptotically approach self-similar ones will exist, and will be 
well approximated by the self-similar solutions that satisfy Eqs.(|17j-pl|) with V{4>) = 0. 



This condition holds for several physically interesting cases. One such case is a massless 
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scalar field in an AdS spacetime background in which we have V = A/8nG. Another is the 
potential for weakly interacting pseudo-Goldston bosons [11]], for which we have 



V(<P) = V cos 2 ^(^y f(p) = 

where V and p are constant, and < p < 1. 

Indeed, as mentioned in the introduction, recent numerical simulations HQ provide 
strong evidence that the critical solution in the case V = A/87iG is continuously self-similar. 
This suggests the curious result that although the cosmological constant is crucial for the 
formation of a BTZ black hole, the actual critical solution does not depend on it. This, 
then, will serve as our motivation to only consider CSS solutions in the following sections. 



Pp- -P) 
4ttG 



III. SOLUTIONS WITH CONTINUOUS SELF-SIMILARITY 

In order to construct CSS solutions, we must do more than assume that Eq.(^) holds in 
the limit as r — > oo. We will assume it as an identity and therefore set V = in Eqs. (|l7| -pl|). 
Doing this, we can drop all dependance on r. We also change variables by choosing to use 
the somewhat more common similarity variable, z = e x . Making these substitutions in the 
field equations, we get 



z 



s" + 2za'(s - zs') = -8irGs(c - zip') 2 , (27) 

s" - 2a' s' = SnGsif' 2 , (28) 

zs" + 2s(za" + a') = 8nGs(p'(c - Zip'), (29) 

s" = 0, (30) 

2zsip" + if'{2zs' + s)- cs' = 0, (31) 

where a prime denotes ordinary differentiation with respect to z. From Eq. (|30|) we find that 

s(z) = a z - b , (32) 

where ao and bo are integration constants. Depending on the values of these two constants, 
the solutions can have very different physical interpretations. Thus, we will consider them 
separately in the following. 



A. ao / and b = 

For the case in which s(z) = aoZ, we find that Eqs. (|27H3lD have the general solution, 

a(z) = -In \z\ + (Jo, ip(z) = cln \z\ + ip l , (b = 0), (33) 

where o"q and ip^ are integration constants with \ = 8tcc 2 G > 0. By rescaling u and v and 
redefining the constant ifl, we can always set ao = 1 and Uq = 0, a condition that we shall 
assume in the following discussion. Introducing the two new coordinates u and v via the 
relations 
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(34) 



we find that the corresponding metric and the massless scalar field are given, respectively, 
by 

ds 2 = -2dudv + v 2/(1+x) d6 2 , 
0=_£_ In (35) 

where is another constant. From Eq. (|35|) we find that 

v -*=(rT3»«- fl - = 8lG ((rf^) 2 «- < 36 > 

which shows that the massless scalar field is equivalent to an ingoing null dust flowing along 
the null hypersurfaces v = const. Because the hypersurface u = is regular, to have a 
geodesically maximal spacetime, we must extend the spacetime beyond this surface. One 
simple analytic extension is to take the range of u in Eq.(|3"5D simply as u G (— oo, +oo). On 
the other hand, the hypersurface v = is singular, and, as shown in Appendix A, the nature 
of the singularity is strong in the sense that the distortion experienced by a freely falling 
observer becomes infinitely large at v = (v — 0). The corresponding Penrose diagram is 
given in Fig. 2. 




FIG. 2. The Penrose diagram for the solutions described by Eq.(|35|) in the text, where the 
double line at v = (i.e. r = 0) represents the spacetime singularity. 



S 



B. ao = and 6q / 



In the case where s( 



-bo, we find that Eqs. 



a(z 



2 



z\ + a. 



have the general solution, 
(ao = 0). (37) 



Again, by a rescaling of u and v, we can set bo = — 1 and ctq = 0. In a manner completely 
analagous to the previous case, we introduce two new coordinates u and v, but defined in 
the current case by 

dv 

du = (1 + y)(— u) x du, dv = r — . (38) 

{1 + X)v x 

Doing so, we find that the corresponding metric and scalar field are given, respectively, by 

ds 2 = -Idudv + u 2/(1+x) d6 2 , 

\n(-u) + (j)l. (39) 



1 + X 



Note that the above metric can be obtained directly from Eq.(p5[) by exchanging the two 
null coordinates u and v. Thus, these two spacetimes must have the same local and global 
properties after such an exchange takes place. In particular, the massless scalar field is again 
energetically equivalent to that of null dust, but this time flowing outwards along the null 
hypersurfaces defined by u = const. Moreover, the hypersurface u = is indeed singular for 
all values of c and, as shown in Appendix A, and in contrast to the claims made in [ 12| , the 



singularity is strong in the sense that the distortion experienced by a freely falling observer 
becomes unbounded as this surface is approached. The corresponding Penrose diagram is 
given in Fig. 3. 



U 




FIG. 3. The Penrose diagram for the solutions described by Eq.(F 
double line at u = represents the spacetime singularity. 



in the text, where the 
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We also note that in |TI| the validity of the solutions of Eq.(^) was restricted to the 
region u < 0, v > 0. However, in order to obtain a geodesically maximal spacetime, as 
we have done here, it is necessary to extend these solutions to the region v < 0. This 
observation will be crucial when we consider the boundary conditions for perturbations in 
section IV. 



C. a b + 



When do^o 7^ 0, it can be shown from Eqs.(|27| [31]) that the constant bo/ao must be 
non-negative, and that the general solution is given by 



a{z) 

(p(z) 



2 X ln 
2c In 



z V4 + az -i/4 



+ 4 



z 1 ' 2 + a 



+ <Po, (a b ^ 0) 



where 



a 



V «0 



(40) 



(41) 



and, as before, x — 8nGc 2 . For a = or a = ±oo, the solutions reduce to those studied in 
the last two subsections. Thus, in the rest of this subsection we will only consider solutions 
with finite, nonzero a. Again, if we rescale the two null coordinates and redefine the constant 
a, we can always set 



a 



-1, C7 1 = -|ln(2 4 ^) 



(42) 



so that at the origin of coordinates, the coordinate u defines proper time. 

In the case that a = 1, the corresponding solutions reduce to those originally found by 
Garfinkle 0. In addition, he found that in the strong field regime the particular solution 
with Xcrit = 7/8 is very similar to the critical solution found numerically by Pretorius and 
Choptuik 0. When a = — 1 the corresponding solutions reduce to those given in [|1^ . 

In terms of u and v, the metric coefficients and the massless scalar field are given by 



a(-u) 1 / 2 + v 



1/2 



r(u, v) = a 2 (—u) - 



-uv) x 



2c In 



a (_ u )i/3 + v i/2 



+ 



-2a 2 e~ 2a < 0, 



from which we find that 

V rV a r 

i.e. V a r is non-spacelike. The curvature is readily found from Eq.(f43D to be 

(— m {\x-W* 

R = R ab g ab = 87rGV a 0V a = 2a X (lQa 2 ) x 



a (_ M )l/2 + u l/2|2(2x+l)' 



(43) 



(44) 



(45) 
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From this, we can see that, for % > 1/2, although the metric coefficients are singular along 
both of the hypersurfaces v = and u = 0, the curvature is perfectly regular there. We 
must therefore extend the solutions beyond these surfaces. 

Note, however that in the case that < x < 1/2, these hypersurfaces at v = and u = 
are singular, and their physical interpretation becomes unclear. For that reason, we will only 
consider the case x — 1/2 in the following. In order to be somewhat systematic in our further 
study of these solutions, we will consider each of the three subcases \ = 1/2, 1/2 < \ < 1 
and x > 1) separately. Finally, it is interesting to note that for all values of x, the spacetime 
is singular at the point (u,v) = (0,0). 

i- x = h 

For x — 1/2) we introduce two new coordinates, u and v, via the relations, 

u = —\/—u, v = \/v, (46) 

and find that in terms of u and v the metric and the massless scalar field are given by 

2 

ds 2 = —-. — r f 2 (u, v)dudv + r 2 (u, v)d6 2 
\a\ 

2/ -\2 -2 

r = a (—u) — v 

= 2cln |/(^)| +0* (47) 



where 



It is straightforward to show 



f(u,v) = a{-u)+v. (48) 



V a rV a r = 8a 2 uve- 2a (49) 



namely, that in Region I where u < and v > 0, the hypersurfaces r(u,v) = const, are 
always spacelike, including the one at r(u,v) = which forms the upper boundary of the 
spacetime [cf. Fig. 4]. 



Using Eq.fl47D to calculate the ingoing and outgoing expansions, we find 



nl \ a \ V nl 2 \ a \ U /rn\ 

^ =-772-' ^ = « 2 ^ (50) 

which are always non-positive in Region I. Thus, all the symmetry spheres 9 = const, are 
trapped spheres (or circles if you prefer) in this region 0. Further, we can show from Eq . (\i7\) 
that 

R = R ab g ab = 87rGV a 0V> = r - ^ l4 , (51) 

[a{— u) + v\ 

which shows that, if a > the spacelike hypersurface r(u, v) = is regular while if a < 
it is singular. From the above expression we can also see that the spacetime is free of 
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singularity on the null hypersurface v = 0, except for the point (u, v) = (0, 0). In order 
to have a geodesically maximal spacetime, we must therefore extend the spacetime beyond 
the surface v = 0. There are, of course, many ways to make such an extension. However, 
we will only consider those that are analytic. By imposing the condition of analyticity, the 
extension becomes unique and is given by simply taking the corresponding solutions to be 
valid in the whole (u, u)-plane. In principle, one might obtain three extended regions in this 
way, I', II, and where /' = {x a : u > 0, v < 0}, II = {x a : u < 0, v < 0}, and 
IV = {x a : u > 0, v > 0} [cf. Fig.l]. 




FIG. 4. The spacetime in the (u, ■u)-plane for the solutions described by Eqs.(|47|) and (^) with 
X = 1/2. The hyper surf aces (p(u,v) = const, are the lines parallel to the one r = in Region I, 
which is singular when a < 0, and singularity- free when a > 0. In constrast, note that the line 
r = in Region II is singular when a > 0, and regular when a < 0. For all values of a, the spheres 
9 = const, are trapped in Region I. In Regions I and II V a is timelike when a < 0, and spacelike 
when a > 0. 



Applying this extension to the above solution, we find that the extended solution in 
Region II is still given by Eqs. (f47D and (|48|), with the corresponding physical quantities 
given by Eqs.(^9|-pl|) but with v < in these expressions. From these equations we can see 
that the outgoing null geodesies are no longer trapped in Region II and that the hypersurfaces 
r(u,v) = const, become timelike. Moreover, from the Ricci curvature it can be seen that 
the surface r(u, v) = in this region is singular for a > while for a < 0, it is regular. In 
both cases, this surface will form the lower boundary of the extended region, such that the 
extended geodesically maximal spacetime in our case will consist of Regions I and II with the 
lines r = as its upper and lower boundaries, as shown in Fig. 4. The remaining portions 
of Fig. 4 are causally disconnected from the extended spacetime. Eq.(|5TD also shows that 
V a is timelike only when a < 0. Thus, in order to interpret the extended solution as 
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representing gravitational collapse of the scalar field, we must take a < 0, a condition which 
we shall assume from now on for the case x = 1/2- The hypersurfaces <f>(u,v) = const, are 
thus straight lines parallel to the singular surface r = in Region I, as shown in Fig. 4. 
The corresponding Penrose diagram is given by Fig. 5. From it, we can interpret Region I 
as the interior of a black hole with the hypersurface v = (or v = 0) as its event horizon. 
In Region II, introducing the coordinates u and v via the relations 

U = — \J — U, V = —y/—V, (52) 

we can find the corresponding metric and massless scalar field in terms of these variables. 
Using a modified definition of z = v/u, we have 



ds 2 = -2e 2a{z) dudv + 
i, v) = cln |— u\ + <p(z) 



2„2, 



-U) S 



z)dd 2 



(53) 



where 



(7 2 



slz 



<f(z) 



In 



^ 1/4 + \a\z 



1/4 



a' — z 
2c In 



a \ + z 



1/2 



+ 



(Region //). 



In the next section, when we study the perturbations of this \ 
form of the metric as the background solution. 

r = 



(54) 

1/2 case, we will use this 




FIG. 5. This is the Penrose diagram for the solutions described by Eqs.(f47|) and ( |48| ) with a < 
or, equivalently, the Penrose diagram for the solutions described by Eqs.(|56|) and (|57|) with a > 0. 
The horizontal line r = is singular, while the vertical line r = is free of spacetime singularity. 
The dashed lines are the hypersurfaces where 4>(u, v) = const. The normal vector to these surfaces, 
V a 4>, is always timelike. In Region I all the surfaces 9 = const are trapped, while in Region II they 
are not. 
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Before moving on to consider the other cases, we note that, if A(u, v) is a solution of the 
Einstein field equations (p^H l5|) with V = 0, then, the solutions 

A'(u, v) = A(ciu,c 2 v), (55) 

also satisfy those equations, where A = {a, r, 0}, and C\ and c 2 are arbitrary constants. As 
an example, let us consider the solution A'(u,v) = A(u, —v), that is, 

2 

ds 2 = —-. — r f 2 (u, v)dudv + r 2 (u, v)d9 2 
\a\ 

It -\2 -2 

r = a {—u) — v 

<j> = 2c\n\f(u,v)\ + (f) 1 (56) 

where 

f(u, v) = a(— u) — v. (57) 

Clearly, this solution can be obtained from that of Eqs. (|47|) and ( 48[) by changing the sign of 
a. Thus, the solution of Eqs.(^) and Q57D with a > also represents gravitational collapse 
of the massless scalar field, and the corresponding Penrose diagram is the same as that given 
by Fig. 5. It can be shown that in terms of u and v defined by Eq. (|52|) the metric and the 
massless scalar field in the present case are also given by Eqs.(p3|) and (0) but now with 
a > 0. 



2. i<x<l 

To begin analyzing Eqs. (^3| -^5|) for the case that 1/2 < % < 1, we again introduce new 
coordinates u and v via the relations 

u = -(-n) 1 / 2 ", v = v x ' 2n (58) 

where we have defined the constant 

n = — - > 1 (59) 
2(1 -X) V } 

(and which will shortly be taken to be an integer, but can be assumed, for the moment, to 
be real). In these coordinates, the metric and the massless scalar field are given by 

ds 2 = — \f(u,v)\ 4x dudv + r 2 (u,v)d6 2 

(16cr) x 

r = a 2 (-u] 2n -v 2n 



where 



2cln\f(u,v)\ + 4>l (60) 



f(u,v) = a(-u) n + v n . (61) 
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The relevant physical quantities are given by 

(16a 2 )* v 2n ~ l 



el 
el 



a 



An r\f(u,v)\ 
2 (16a 2 )* {-u) 2n - 1 



An r\f(u,v)\ 



V„rV°r = -(16a 2 Y ' '"' ) "~ > 



-uv) n 1 



a(-u) n + v n \ 4x 

R = 8vrGV a 0V> = 2a X (16a 2 )* \ ^ 7 , } . (62) 

Because we have n > 1, we can see from the above expressions that all the surfaces of 
symmetry, = const, are trapped in Region I and that the spacetime is bounded from 
above by the hypersurface v = |a| 1 / n (— u) {i.e. r = 0), which is singular for a < 0, and 
regular for a > 0. Note that as v — > + , we have R — > and the spacetime is perfectly 
regular on v = as well. Therefore, in order to have a geodesically maximal spacetime, 
we must extend the spacetime beyond this surface to v < 0. As in the previous case with 
X = 1/2, we can simply take the above solutions to also be valid in Region II. However, 
this extension will not be analytic unless n is an integer. When n is not an integer, the 
extended metric will not even be real. As an example, consider the case where n is rational, 
i.e. n = (2m + 1)/(2Z), where I and m are integers. A possible extension in this case would 
be 

A n (u,v) = A\u,-v). (63) 

Clearly, this extension is not analytic, too, but it guarantees that the metric in Region II is 
real and the extension has the maximal order of derivatives in comparing with any of other 
extensions. 

However, as mentioned previously, we will only consider those cases where the extensions 
are analytic. Therefore, we shall restrict ourselves in the following to those cases where n is 
an integer. For such an extension, the relevant physical quantities in Region II will also be 
given by Eq . (|62|) but now with v < 0. With this, we find that in Region II, in order to have 
V a be timelike we must have 

U \2Z + 1, «<0, [m) 

where I is an integer. To further study these maximally extended solutions, let us consider 
the cases a > and a < 0, separately. 

Case (a) a > 0: In this case we must choose n = 21 so that V a is timelike in Region II, 
and the corresponding solutions can be interpreted as representing gravitational collapse in 
this region. V a is spacelike in Region I, while on the hypersurface v — it becomes null, as 
we can see from Eq . (|62|) . The hypersurfaces <p(u,v) = const, are shown in Fig. 6. Similar 
to the last case, we may interpret Region I as the interior of a black hole, since the surfaces 
9 = const, are trapped, although no spacetime singularity is present in this case [except for 
the point (u, v) = (0, 0)]. This black hole can be considered as the final state of the collapse 
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of the massless scalar field in Region II. This can be seen most clearly from an analysis of 
the energy-momentum tensor. Near the null hypersurface, v = 0, the only non-vanishing 
component of the energy-momentum tensor is given by 



T- 

-L it 



(65) 



which represents a pure energy flow along the null hypersurface u = const, moving from 
Region II into Region I. The corresponding Penrose diagram is also given by Fig. 5, but 
now both of the vertical and horizontal lines r = are free of spacetime singularities, except 
for the point (w, v) = (0, 0). 




FIG. 6. The hypersurfaces 4>(u,v) = const, in the (u, v)-plane for the spacetime described by 
Eqs.(|0| HIOHII) with a > and n = 2/ (q < and n = 2/ + 1), where I is an integer. The 
normal vector to these surfaces X7 a (j) is timelike in Region II and null on the hypersurface v = for 
both of the two cases a > and a < 0, while in Region I it is spacelike for a > and timelike for 
a < 0. For both cases, the spheres 6 = const, are trapped in Region I but not in Region II. The 
spacetime is free of singularities on the surfaces r = in both regions when a > 0, while for a < 0, 
the surface r = in Region I is singular. The sole exception to this is the point (it, v) = (0, 0) 
which is always singular. 



In Region II, introducing u and v via the relations, 

= — (— v) Al , (a > 0, region //) 



u 



(66) 



we find that the metric and the massless scalar field in terms of u and v (again using the 
modified definition of z = v/u) are given by 
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ds 2 = -2e 2a{z) dudv+(-u) 2 s 2 (z)d0 2 
(u, v) = cm | — u\ + ip(z), 



(67) 



where 



2 V4 + | a U-l/4 



s(z) = a — z 



(f(z) = 2c In 



\a \ + z 



1/2 



!>o, (region II) . 



(68) 



In the next section, when we study perturbations we shall use this form for the metric as 
the background spacetime for this case. 

Case (b) a < 0: In this case we must choose n = 21 + 1 so that V a is timelike in Region II. 
This region is bounded by the timelike hypersurface r = 0, on which the spacetime is regular. 
Unlike the case a > 0, in Region I now V a is still timelike, while on the hypersurface v = 
it is null. The hypersurfaces 4>{u, v) = const, are also given in Fig. 6. However, now the 
spacelike hypersurface r = in Region I is singular and the surfaces 9 = const, are trapped 
surfaces. Near the null hypersurface, v — 0, the energy- momentum tensor also has only 
one non- vanishing component, given exactly by Eq.(^). Therefore, in the present case the 
corresponding solutions can also be interpreted as representing gravitational collapse of a 
massless scalar field, with a black hole finally formed in Region I. The corresponding Penrose 
diagram is again that of Fig. 5. 

In Region II, on defining u and v as 



u 



-(-u 



>2(2/+l) 



-l-V 



,2(2Z+1) 



(a < 0, region //) 



(69) 



we find that in terms of u and v the metric and the massless scalar field are exactly given 
by Eqs.(|67|) and (|68|). 



3. X>1 

In this case it can be shown that the spacetime is already geodesically maximal in the 



region v > 0, u < 0, and does not need to be extended beyond the hypersurface v = |L2 
Indeed, it is found that the null geodesies u = const, have the integral 

near the hypersurface v — 0, where rj denotes the affine parameter along the null geodesies, 
and f3 and r/o are integration constants. Thus, as v — > + , we always have 7] — > ±oo for 
X > 1. On the other hand, from Eq.(^3) we can see that the spacetime is singular on the 
spacelike hypersurface v = —au (or r = 0) for a < 0, and regular on this surface for a > 0. 
Moreover, within the entire region, the symmetry spheres 6 = const, are trapped surfaces, 
as can be seen from the expansions 6i and 9 n which are always non-positive 

1 2 

6i = --e~ 2a , e n = -ye~ 2 °, (u<0,v> 0). (71) 
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It is therefore difficult to interpret the solutions in the present case as representing gravita- 
tional collapse. 




FIG. 7. The Penrose diagram for the solutions given by Eqs.(|74|). The dashed lines represent 
the hypersurfaces <p{u,v) = const., the normal vector of which is always timelike. The whole 
spacetime is free of singularity except for the point (u,v) = (0,0). 



However, applying Eq. ( p3|) to the above solutions, we can consider the solutions in Region 
II, but which are now causally disconnected from Region I. These solutions are given by 



a(-u) 1/2 + (-v 



1/2 



4x 



(16a 2 )* 
r(u, v) = a 2 (—u) + v, 



[uv * 



(u, v) = 2c In 



a 



+ { _ v y/2 ( U ,V<0). 



From these expressions, we find that 

0, = - a - 2 e„ 



-2a 



(72) 



(73) 



which is non-negative; that is, the outgoing null geodesies are now no longer trapped. One 
can show that the null hypersurface v — is also future null infinity of Region II. Thus, the 
spacetime is geodesically maximal in this region. As a result, we also have 
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[uv 



lX-1/2 



R = BirGVaW* = -2a X (16a 2 )* - ^' 2(2x+1) , (74) 

which shows that, when a < the spacetime in this region is singular at r = 0, and when 
a > the spacetime is free of such singularities. Moreover, from Eq.(|74D we can also see 
that the normal vector to the hypersurfaces (f)(u,v) = const, is always timelike when a > 0, 
and quite similar to those given by Fig. 6. Thus, the solutions with a > can be thought 
of as representing gravitational collapse starting from past null infinity at u — — oo. Thus, 
we will only consider solutions with a > for the current case of x > 1. As a result, 
the spacetime is free of singularities at the beginning of the collapse, but with a spacetime 
singularity forming at the point (u, v) = (0, 0). The corresponding Penrose diagram is given 
in Fig. 7. It is interesting to note that in this case, although no event horizon is formed, the 
spacetime singularity at the point (u, v) = (0, 0) is not naked and the only observers who 
are able to see the singularity there are those observers who arrive at that point. 

It can be shown that in this case (with a > 0) the metric and the massless scalar field 
in Region II can be also written in the same form as that given by Eqs. (|57D and (|B5[) . 

In fact, from Eqs. ([53|-|54]) and Eqs. (|67H68|) , we find that in the case a &o 7^ 0, all the solu- 
tions that represent gravitational collapse of the massless scalar field can be cast in the form 
of Eqs.fflTfy and §68J in Region II. For that reason, when we consider linear perturbations in 
the next section, we will always refer to Eqs. fl67|) and (|68| ) with x > 1/2 when considering 
the case of ao&o 7^ 0. In addition, because it is only the absolute value of a which enters 
these equations, we will from now on and without loss of generality assume that a > 0. 



IV. LINEAR PERTURBATIONS OF THE SELF-SIMILAR SOLUTIONS 

We now turn to a study of the linear perturbations of the self-similar solutions presented 
in the last section. Let us consider perturbations given by 

a(z, t) = oq(z) + ea 1 (z)e kT , 
s(z,t) = s (z) + esi(2)e fcr , 

ip(z } r) = ip (z) + tip l (z)e kT } (75) 

where e is a very small real constant, and quantities with subscripts "1" denote perturbations, 
and those with "0" denote the self-similar solutions, given, respectively, by Eq.(|33D with 
s = z , al = 0, Eq.([T?p with s = 1, a* = 0, and Eqs. (§70 and (§§ with x > V 2 - It 
is understood that there may be many perturbation modes for different values (possibly 
complex) of the constant k. The general perturbation will be the sum of these individual 
modes. Those modes with Re(k) > grow as r — ► oo and are referred to as unstable modes, 
while the ones with Re(k) < decay and are referred to as stable modes. By definition, 
critical solutions will have one and only one unstable mode. It should be noted that in 
writing Eq. ([Top , we have already used some of the gauge freedom available to us to write 
the perturbations such that they preserve the form of the metric (0). However, this does not 
completely fix the gauge freedom and we will return to this issue when we consider gauge 
modes. 

To first order in e, it can be shown that the non- vanishing components of the Ricci tensor 
are given by 
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R, 



If 

y,uu 

r 

„2t 



uIsq 



z 2 s' ' - 2za' (zs'q - s ) 

{2+k)r ( 

z 2 s'[ + 2z{k - za' )s' 1 
- 'z 2 So 1 (sg - 2a' s' ) + k{2za' + l-k) 
-2(zs' - s Q )(za' 1 + ka x ) \ 



ufco 



Re 



,2r 



U^Sq 



; (2 + fc)T 



(r uv + 2ra, 



•2t 



UqSq 



zs'q + 2s (za' ' + a' 

e (2+k)r < 



tt§s 



zs'l + 2zs a" + ks[ + 2(1 + k)s a[ - zSq Sq'si 



2e~ 2<T rr 



2zs s e 



uv 
//_— 2<jo 



(76) 



To zeroth order the Einstein field equations are given by Eqs.(|T7]-|20|), while to first order 
they take the form 



Si 



z 2 s" + 2z{k - za' )s[ 

- ~z 2 s 1 {4 - 2a' s' ) + k{2za' + l — k) 

- 2(zs' - s )(za[ + kax) 

s" - 2a' s[ - s 1 (s ' - 2a' s' ) sx - 2s' a[ 
zs'l + ^zs a'( + ks[ 

+ 2(1 + k)s al - zsq 1 s'qS 1 
zs'{ + k s 'i + zs o ( s o * s i ~~ 2<xi 



-167rGs (z(p' - c) (zip[ + kipx) (77) 

— lGnGsoipQipl (78) 

-8ttGs l(2z<fi'o ~ c) ip[ + ^oV^i] (79) 

0. (80) 



To zeroth order, the equation for the massless scalar field is given by Eq.(gjj) while, to first 
order, the scalar field equation is found to be 



2zs (p'l + [2zs' + (1 + 2k) s ] f'x + ks' ifx 

+ (2zcp' - c) s[ + [2z(f% + (1 + k)<p' Q ] Sl = 0. 



(81) 
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In order to be somewhat systematic in our examination of the perturbation equations, 
we will consider separately the same three cases that we did in section III, namely a® 7^ 
0, b = 0; a = 0, b 7^ and a b 7^ 0. The notation will of necessity become rather 
complicated and we ask for the reader's indulgence. For our part, we will try to keep the 
complexity under control and the physical interpretation near the fore. 



A. ao / and bo = 

In this case, the self-similar background solutions are given by Eq . (|3~3"D with s (z) = z 
and o\ = 0. Substituting these solutions into Eqs.([f7|-[81|) we find 

zs'[ + ks[ = (82) 

(k - X ) [zs'i + (k — = (83) 

z{k + xWi - XSi + 2z 2 a' 1 = 16ncG z 2 ip\ (84) 

z 2 a" + (1 + k)za[ = -AixcG [ktp x + z<p[) (85) 

2z 2 <p" + (3 + 2k)z(p[ + k(pi = -- [zs[ + (k - l) Sl ] . (86) 



It should be noted that in writing Eqs.(|83|-pq), we have used Eq. (|82|) . As usual, Eqs.(p3|) 
and (|J) are constraints, and Eqs. (|82|) , (jS5j) and (|^) are the dynamical equations. Due to 
the redundancy in the Einstein equations, these are not all independent and, e.g. it can be 
shown that Eq.(^) can be written in terms of the others. 

Once these equations are integrated, one needs boundary conditions to complete the 
solution. Such boundary conditions are crucial for a proper determination of the eigenmodes 
but are rather subtle in general relativity. Indeed, there is no fixed prescription to follow 
|7|JT^|. In the present case, the background spacetime has four boundaries, u = ±00 and 



v = 0, 00. However, since the perturbation equations are of second order, it is usually 
sufficient to impose boundary conditions on only two hyper surf aces. As a x , s\ and ifi are 
functions of z only, it is rather natural to choose the hypersurface z = Zq as one of the two 
boundaries, that is, we assume that the perturbations are turned on at z = z , where z 
is an arbitrary constant. Prior to z = Zq the perturbations are zero. In order to match 
the perturbations smoothly across this surface, we require that S\(zq) = = ^1(20). Also, 
by using the coordinate tranformations (D) we can always set 01(2:0) = 0. Then, on the 
hypersurface z = z we have the boundary conditions, 

<7i(zo) = si(z ) = cpi{z ) = 0. (87) 

The second boundary has to be one of the two hypersurfaces, v — and u = 00. As 
shown in the last section, the hypersurface v — is always singular in the present case 
and we do not know how to impose boundary conditions on such a surface.f] Therefore, in 



4 Sometimes the condition is imposed that the perturbations should be "less singular" than the 
background. However, it is not always clear if the expression "less singular" should refer to the 
spacetime curvature or to the metric components. 
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the following we shall choose the hypersurface u = oo as the second boundary. Because 
this surface represents future null infinity for this spacetime, it is reasonable to assume 
that the perturbations grow slower than the background. In the present case, we have 
s (z) ~ z, a (z) ~ ln(z), (fio(z) ~ ln(z). However, considering the fact that e r ~ (— it) -1 , 
we might require that, as u — > oo, 5r(u,v) = e(—u)s\(z)e kT and 5ip(u, v) = e{pi(z)e kr be 
finite while 5a(u, v) = eai(z)e kr is requied to grow slower than ln(— u), as u — ► oo. Instead, 
we will choose to impose somewhat stronger conditions here, namely 

Sr(u, v), 8(p(u,v), Sa(u,v) — > finite, (88) 

as u — > oo. Clearly, if we find unstable modes with this stricter criterion of regularity, then 
in general the solutions will still be unstable. 

Before proceeding to the full perturbation analysis, we need to address questions con- 
cerning gauge modes. Under the coordinate transformations 

u— >w + e^i(w), v — > v + e£ 2 (v), (89) 

we generate perturbations of the form 

^ = dx +(> 

2 [ [ f u J 

5r = e£ 2 (u) 

^ = fC f§M_«iM), (90) 

where ^i(it) and 6( w ) are arbitrary functions. For such perturbations to take the forms of 
Eq. fl75|) , we must choose 

= ^VH*) 1 "*, = (9i) 

where c\ and C2 are arbitrary constants. Then, the solutions to the perturbation equations 
due to the gauge transformations are given by 

si(z) = c 2 z 1 - k 

a l (z) = C j(l-k + X )z- k -^ c (l-k-x) 

ipx(z) = ci + cc 2 z~ k . (92) 

It can now be seen that our boundary conditions, Eqs. fl8"T|) and (|58] ) eliminate each of these 
gauge modes. 

Now let us turn to solving the perturbation equations given by Eqs.(|52}-|8"6|). From Eq.(|53D 
we can see that there are two subcases which we need to consider separately, k = x an d 
k ^ x- 
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In this subcase, we have 

k = x = 87rGc 2 > 0, 
that is, the constant k is real and strictly positive. Integrating Eq.( 



we obtain 



«i0) 



/31n|z|+s?, k 



X 



(3z 



i-x 



+ 4, 



(93) 



(94) 



where (3 and s® are real integration constants. 



Case (a) k — x — 1 ; From Eqs.(^-|SB|) we find that the solutions to the perturbation 
equations can be written 



si(z) = (3ln \z\ + s[ 
o~i(z) = 87icG <fi(z) 



C!Z~ 1/2 + c 2 z~ x + — In \z\ 

z 



(95) 



where ci, C2 and 0"° are additional real, integration constants. Applying the boundary 
conditions at z = z , Eq. fl8~T|) , to these solutions, we find that 

5° = — j3\nz 



i = — 

Cl : 



2z 

1/2 



-2, 







c 2 + c/?ln(^o) 



(96) 



from which we obtain 



r(u, v) = (— u) sq(z) + esi(z)e 



v + e(3uQ In 



kr 



-ZqU 



(97) 



From this form for r(u,v), we can see that the condition of finiteness as u —>■ oo, Eq.(^), 
leads to (3 = 0. Setting (3 = in these solutions yields that 5a and 5ip are all finite as 
u — > oo. Therefore, in this case (A; = %), there is a perturbation mode with k = 1 which 
satisfies the appropriate boundary conditions, Eqs.flSTD and (^). Of course, because k > 0, 
this is an unstable mode. 



Case (b) k — x 1 ; I n this case, the perturbation equations given by Eqs. (]82|-j36D can 
again be integrated to yield 



sAz 



f3z l - x + s° 



<r x (z) = SncGMz) ~ f (2x - 1) ^ X " 4^ + 



-1/2 



„-l 



+ C2Z x + CSV Z 



CiZ l l 2 + C 2 Z 1//2 In |z| + CsJ 2! 



„-l 



Ml/2, 
k — 1/2. 



(98) 
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On application of the boundary conditions, Eqs. (|87| ) and fl38|), at z = z and u — > oo, it can 
be shown that there are nontrivial perturbations for k > 1/2 while for k < 1/2, no solutions 
exist which are compatible with the boundary conditions. 

We conclude this section by noting that all the solutions given by Eq. (|33"D of section IIIA 
with x > 1/2 have an unstable mode with k = x- 



2. k + X 



In this subcase, the constant k can be complex, and from Eqs. (|82|-|36"D we find that we 
can again integrate the perturbation equations yielding 

s 1 (z)=(3z 1 - k 

<n(z) = 8ttcG Vi (z) + ^/3(1 — k — x)z~ k + <7? 

{z) _{c lZ -^ + c 2 z-\ Ml/2 { ) 

where (3, a®, Ci, and c 2 are additional arbitrary (but now possibly complex) constants of 
integration. It is interesting to note that setting c\ = = s® and c 2 = cj3 the above 
perturbations will reproduce the gauge modes given by Eq.(|9"2"D that do not satisfy the 
boundary conditions (p7[) and d88|) . However, the general perturbations represented by 
Eqs.([99|) will satisfy these conditions, provided that 

P = a° x = 

Ci — — C 2 2 

i?e(A;)>^, (A;^ X ). (100) 



Because Re(k) > 1/2 > for any given x, we can see that the solutions given by Eq . (p3|) 
are unstable with an infinite number of unstable modes. This result is consistent with 
King's conjecture |13| , namely: Non-scalar singularities are unstable to perturbations and 



will eventually give rise to scalar singularities. 



B. ao = and b$ ^ 

In this case, the self-similar background solutions are given by Eq. (|37D with sq = 1, a-g = 
0. Substituting these solutions into Eqs.([77|-|8TD we find that 



z 2 s'[ + {2k + x)zs[ - k(l — k - x)si (101) 

+ 2(^0"! + za[) = WncGikipx + z(p[), (102) 

za'l + (1 + k)a[ = AncGip'^ (103) 

zs'l + X s[ = 0, (104) 

zs'[ + ks\ = 0, (105) 

2z<p" + (1 + 2k)ip[ = cs[. (106) 
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Note that in writing the above equations, we have used Eqs.( [LQ4] ) and ( |105| ) to simplify the 
expressions somewhat. 

Following arguments similar to those given in the previous subsection for ingoing solu- 
tions with a ^ and bo = 0, we here impose boundary conditions on the hypersurfaces 
z = zq and v = oo. The conditions on the hypersurface z = zq are the same as those 
given by Eq.(^). On the hypersurface, v — oo [z — oo), we take as our condition that 
the perturbations should grow more slowly than the unperturbed, background solution. At 
future null infinity these conditions become 

Si(z), ifi(z) — > finite, 

&i{z) grows slower than lnz, (107) 



as z — > oo. As a result, it can be shown that the gauge modes of Eqs.(|90 H9ll) are given, in 
the current case, by 



ai(z) = c z k + ci, s x (z) = ipx(z) = 0. (10£ 



Returning to the perturbation equations in Eqs. ( pL0lHl06|) , we can see from these that 



there are two distinct possibilities, s[(z) ^ and s[(z) = 0. In the former case it can 
be shown that the above equations have no solution consistent with the above boundary 
conditions at z = z and v — ► oo. Thus, in the following we need only consider the case 
s[(z) = 0. In this latter case, the perturbation equations can be rewritten as 

si(z) = si, (109) 

za \ + ka 1 = %ncG{k Vl + zip[) + -s?/c(l - k - x), (110) 

2z<f/-l + (l + 2k)<p' 1 = 0. (Ill) 

The general solution to these equations is then 

*i(*) = a? 

a^z) = tocGniz) + a° lZ - k + ^(l-k- X ) 

^ z) = \(3z^y* + rt k^ l ^ 

where /3, ip\ and a® are integration constants (real for k = 1/2 and possibly complex other- 
wise) . 

Note that the gauge modes of Eq. (|108|) can be obtained from Eqs. (|112|) by setting a® = c 
and si = (3 = 0. However, these gauge modes are completely eliminated by our boundary 
conditions (0) and fllDTD . 

Applying the boundary conditions to the general perturbation solutions given by 
Eqs.( |112| ) eliminates all the perturbations for k — 1/2 while for k ^ 1/2 there are non- 
trivial perturbations which do satisfy the boundary conditions provided we have 

s° = <x° = 

<p° = ~[3zt 2k)/2 
Re(k) > ~. (113) 
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Because the nontrivial perturbations in this case have Re(k) > 1/2 > 0, we can see that 
all the solutions in section IIIB given by Eq. (|37D are unstable with an infinite number of 
unstable modes. 



C. a b + 

With this most general case, we can interpret the solution as collapsing and might hope 
that by examining the perturbative structure, we can determine if there are critical solu- 
tions. The self-similar background solutions are given by Eqs.(|67|) and (68) with \ > 1/2- 
Substituting these solutions into Eq.(BO) we find that 



zs'[ + ks[ = 0, (114) 



which has the general solution 



where (3 and s\ are the integration constants. 

Before solving the rest of the differential equations, fl77j-fT^) and (|STD, let us again first 
consider the boundary conditions. As we can see from Figs. 5 and 7, Region II of the 
spactime has only three boundaries: the center of the spacetime at z = a 2 , the event 
horizon at z — 0, and past null infinity at u — — oo. In this paper we choose z = a 2 and 
z = as the surfaces on which boundary conditions will be imposed. 

To see what kind of boundary conditions should be imposed at each of these surfaces, 
first notice that z = a 2 is the center of the background. It is therefore natual to assume 
that this will remain true even after the spacetime is perturbed, that is, s\(z = a 2 ) = 0. In 
addition, because the background is free of spacetime singularities at the origin in Region II, 
the perturbed spacetime should remain regular there as well. It can be shown that this will 
only be the case if (T\{z) and the quantity [(a 2 — z)(p[ — 2ktpi] remain finite. Summarizing 
the above, we impose the following boundary conditions at the origin of coordinates 




Ul-y)^^--2k<p 1 (y)} ~ finite, (116) 

where we have defined the new variable y = z/a 2 such that z = a 2 corresponds to y — 1. 

The other boundary at z = 0, or v = 0, represents the event horizon.^ We would like 
to be able to assume that no matter comes out of Region I into Region II. Naively, we can 
express this as 



5 Note that for the case \ > 1, the hypersurface v = is marginally trapped, as one can see from 
Eq.(|73|). However, in our discusssion of boundary conditions, we shall treat it as an event horizon 
too. 
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<t>,v = 



(117) 



as v — > . However, it can be shown that this condition does not hold even for the 
background solutions. Indeed, as v — > 0, we have 

dcf>o(u,v) c , . 

But this is nothing more than the statement that in the coordinates u and v , the metric 
coefficients are singular there. Having found the analytic extension of these collapsing so- 
lutions across the hypersurface v = in the coodinates u and v, we should understand the 
no outflow condition, Eq.( |117|) , in terms of v rather than the "bad" coordinate v. This 
argument is further justified by noting that for the unperturbed, background solution, we 
have (po,v — > 0, as v — > 0. In terms of our rescaled self-similar variable y, it can be shown that 
this condition is equivalent to y x ip' x {y) — > 0. This then, becomes our boundary condition on 
<p(y) at y = 0. 

For the boundary conditions on the other fields at y = 0, recall that when we extended the 
background solutions from Region I to Region II, we required that the extension be analytic. 
Otherwise, the extension will not be unique and in some cases may even be meaningless. To 
have consistency between the unperturbed and perturbed spacetimes, we will assume that 
the perturbations are also analytic across y = 0. Given this, we find that the boundary 
conditions at the event horizon y = become 

si(y), <ri(y), <pi(y) are analytic, 

y*^Uo. (119) 
dy 

Again, a word about gauge modes is in order. It can be shown that in the present case 
the gauge transformations (|89"D lead to the perturbations 

i-fc 



.,(») = ,3 (a 2 y) ~ +s° 



Let us now consider the remaining perturbation equations given by Eqs.(|77H79|) and (pT|). 
It turns out to be convenient to consider the two cases k — 1 and k ^ 1 separately. 

1. k = 1 

In this case from Eq. (|115|) we find that 

s l {y)=(3\n{y). (121) 
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where, in writing the above expression, we have chosen = —2/3 ln(a) so that the boundary 



condition si(y = 1) = is satisfied. Inserting Eq. (|121[) into Eqs.(|77|-ff9|) and (pTD , we find 
that only two of the four equations are independent 



a 1 (y) = 8ircG (l - y 1/2 
P 



y 1/2 [i + y 



1/2 



d<Pi{y) 

dy 



+ <pi 



2a 2 y{l + y 1 / 2 ) 2 



(i-x)(i + y 2 ) 
+ 2(i + x )y + '2y 1/2 (i + y) 



(122) 
(123) 



where 



1 a 1 x 3 

a = l, 6=2, A = -, 



/(y) 



c/3 



2a 2 2 /(l + y 1 /2)2 



(124) 



Eq. ( |123j ) is the inhomogeneous hypergeometric equation 14| , and the corresponding general 
solution is given by 



My) = ¥i(y) + v s M, 



(125) 



where <p\(y) is a particular solution of the inhomogeneous equation ( |123|) , and <Pi(y) is the 
general solution of the associated homogeneous equation, 



y(l-y)^ + [X-(a + b+ l)y] ^ - abcp, = 0. 



(126) 



The general solution <Pi(y) can be written as a linear combination of two independent solu- 
tions <p\ and ip\ 



<Pi{y) = ci <Pi(v) +C2<Pi(v), 
where c\ and c 2 are two arbitrary constants, and (p\ and <p\ are given by 



(127) 



2'2- y 
1 1 
•2' 0i 2 



v 2 1 = y- 1 / 2 F(-,0;-;y) = y 



1/2 yi_ y i/2 

-1/2 



(128) 



with F(a, b; A; y) denoting the hypergeometric function. A particular solution of the inho- 
mogeneous equation ( |123| ) can be found and is given by 



c/3 



a 2 {l + y x l 2 ) 



y 



1/2 



(129) 
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From Eqs.( |124] - 129 ), we can show that as y — > 1, we have 



y 



1/2 



Ci In 



1 + y 



1/2' 



y 



1/2 



+ c 2 



(130) 



As a result, the boundary condition for ipi in Eq.( 116 ) holds only if c\ = while the condition 
for ai(y) at y = 1 does not impose any additional restrictions. 

With regard to the boundary conditions at y = 0, we find from Eqs.( |12l1 ) and (|129|) that 
the conditions on si(y) and <7i(y) further demand c 2 = (3 = 0. Therefore, in the case that 
k — 1, the boundary conditions ( |116| ) and (|119|) eliminate all the perturbations. 



2. k + 1 



When fc ^ 1, we have 



si(y) = c ct 2 (y 



i-fc 



1 



;i3i) 



where Co = a~ 2k /3. Note that now the constant k can be complex, as can cq. Substituting 
Eq.( |131| ) into Eqs. (|77| -f7^) and (^T|), we find that now there are also only two independent 
equations, which can be written as 



kai(y) 



X, 



i-y 



1/2^ 



,1/2/ 



1/2 ^ , , 



co(i - fe) 

2(1 + ^/2) 



2/C 1 - 2/)^^" + [A - (« -h 6 -h 1)2/] 



(* + X ) (i + y 1/2 ~ fc 
+ (* - x) (y 1/2 + 

dip 



but now with 

a = k 

m 



dy 2 



»4 



dy 



aVi = /(!/)> 



(132) 
(133) 



c/3 



2a 2fc |/ fc + 1 /2(i + y V2)2 



(l-^l-j/V 2 "^-/" 1 ) 



(134) 



Similar to the previous case, the general solution of Eq.( |133| ) can be written in the same 
forms as those given by Eqs.( |125| - |127|) , but <p\ and (p\ now are given by ip\ = F(^,k;k+^;y), 
and iff = F(^, k\ 1; 1 — y). In addition, it can be shown that a particular solution of the 
inhomogenuous equation ( |133j ) is given by 



c(3 l + y l ' 2 - k 



a 2k i +J/ i/2 ' 

Thus, the general solution for <fi(y) in the present case is given by 



(135) 
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<pi{y) = d F(i + + c 2 F(l fc; 1; 1 - y) - ^ V/iV " (136) 
2 2 2 a zlt 1 + y 1/z 

From Eqs.( p.31| ), ( |132|) and ( |136| ) we can see that when c\ = c 2 = 0, the perturbations reduce 
to those given in Eq.( |120|) and which are purely gauge modes. Therefore, in the following 
we will restrict ourselves to those solutions for which |ci| 2 + | C2 1 2 7^ 0. 

Considering first the boundary conditions at the origin, y — 1, we note the following 
limiting relations as y — *■ 1 

F(i k; k + i; y) - -±(2fc - l)A(fc) ln(l - y) + A a (y; fc) 
F(|, A; + 1; * + ?; y) -> ^(4fc 2 - l)A(fc) (^7^) + ln (l " *)) + Mv\ *) ( 137 ) 



where A(k) = r(& — |)/ (T(A;)r(~)J, and Ai(y; fc) and A 2 (y; k) are the finite pieces of the 
hypergeometric functions in the limit as y — > 1. Then, from Eqs.( |136| - |137| ) and the relation 



^-F(a, b; A; 2) = ^F(a + 1, 6 + 1; A + 1; z), (138) 
02 A 



we find that the boundary conditions given by Eq. ( |116|) at y = 1 will hold, provided that 



ci = 0. (139) 

For the boundary conditions at y = as given in Eq.( |119| ), we first recall that our 
imposition of analyticity on the extended spacetime exploited the barred coordinate system. 
Doin the same here, we notice that y = a~ 2 v/u = a~ 2 (v/u) 2n and 

Sl {y) = c a 2 a" 2 U -1 . (140) 



Thus, to have si(y) be analytic across the boundary y = (i.e. v = 0), we must have 

k = — <l, (141) 
In 

where m is a positive integer.^ Due to the properties of the hypergeometric functions, it is 
convenient to consider separately the cases k = I + 1 and k 7^ 1 + 1, where Z is a non-negative 
integer. 

Case (a) k — I + ^: In this case, as y — > 0, we find that 

( 2' fc ' i ' i yj r(i + or(|)^ (i-Opp! y 

j-lny. (142) 



r(I-0r(|) 



6 In principle, m could be a negative integer as well, but our interest here is in unstable modes 
with k > 0. Hence we restrict to m > 0. 
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Inserting the above expression into Eq. (|136|) and considering Eq. ( |139|) , one quickly sees that 
the condition that <fi{y) be analytic across the hypersurface y = can only be true for 
k < 1/2 and c 2 = 0. Because the boundary conditions as y = force us to take C\ = c 2 = 0, 
we can see that the corresponding perturbations are purely gauge. Thus, we need no longer 
consider the case of k — I + ^. 

Case (b) k I + h In this case, it can be shown that in the limit y — > 



F{~, k + 1; 2; 1 - y) -> k~ l A{k) (2* - l)y~ l/2 ~ k - (1 - % 1/2 " fc 



+ A 3 (y;A;) (143) 



where A 3 (y; k) is again the finite portion of the hypergeometric portion in the limit as y —>■ 0. 
Inserting this into the general solution for <^i(y) as given in Eq. ( |136| ) together with the fact 
that we still have c\ = 0, we can find the limiting behavior of the derivative of phi\(y) as 
y — > 0, namely 



dy 



{cc 



-1/2 



+ 2(1-%" 



+c 2 {2k-l)A{k)y- 1 ' 2 - k 

+ [cc (10k - 7) - 02(1 - k)A(k)] y l ' 2 - k ) 

+A A (k) 



(144) 



where, again, A±{y; k) is the finite part. Now, one can see that the boundary condition of 
vanishing of y x <£[ (y) at y = can only be true if 



k < y and y > - 

2 2 



(145) 



because we want to demand that c 2 ^ 0. Otherwise, the perturbations would again be 
purely gauge. 

We are not yet done as we must still consider the effect of imposing analyticity on o\ (y) 
at y = 0. From Eq. (|143|) and the asymptotic behavior of 



F(-,k;l;l-y)^A(k)y 1 / 2 - k + A 5 (k), 



(146) 



X i A n_\ _ ^ _ . 1/2— A; , \ 



<*i(y) -> - ( A ( k )°2 - cc ) y 
c 



,-k 



err 



(147) 



as y — > 0, we get that 

[A(%(2A:-1)c 2 -(1-x-A:)co" 
as y — > 0. Thus, to have o"i(y) be analytic across the hypersurface y = 0, we must assume 

(148) 



k < - and c 2 = — — , , , en. 

2 X (2A;-1)A(A; 



When k = \ — \ [i.e. m = 1), we can see from the above expression that c 2 = 0. In other 
words, the mode m = 1 is a pure gauge mode. Thus, physically relevant perturbations will 
have m/1. 
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Combining Eqs. (|61|), ( |141|) and ( pL45| ) with Eq. (|148|) , we find that in the present case all 
the perturbations with 

n(l-m)r(fc)r(i) 
01 ' ° 2 ( m -n)(2n-l)T(k-l) C0, 

Tfl 

k = — , l<m<n— 1, n > 1, (149) 

satisfy the boundary conditions ( |116| ) and ( |119| ). 

Therefore, for any parameter n > 2, the corresponding solution of Eqs. fi(F\ ) and ( ^6Bj ) has 
(n — 3) unstable modes. In particular, the one with n = 4 only has one unstable mode. By 
definition this solution represents a critical solution ||. This is consistent with Garfinkle's 
original observations @. However, it is a departure from the perturbation analysis of || 
where they calculate the n = 2 solution to have a single unstable mode. The discrepancy 
between that work and ours would seem to lie in the choice of boundary conditions for 
the perturbation problem. Though both works consider the same boundaries at the origin 
(y = 1) and the horizon (y = 0), we have the added condition that nothing escape from the 
horizon. This would seem to eliminate some of the additional modes found in PI. 



V. CRITICAL COLLAPSE OF SCALAR FIELD WITH POTENTIAL 



In this section, following |L0] we argue that the critical solution with n = 4 (or x = 7/8) 
for the massless scalar field collapse found in the last section is also the critical solution for 
the collapse of a scalar field with any potential V((f>) and cosmological constant A, provided 
that the condition (|26| ) holds. As shown in section III, because of the presence of these terms, 



globally self-similar solutions do not exist, although in the limit r — > oo asymptotically self- 
similar solutions do exist. Fortunately for the study of critical collapse, this is sufficient. 
Indeed, initial configurations for collapsing gravitational systems can be quite different. This 
is the case even in situations where globally self-similar solutions are permitted. Thus, in 
generic cases the spacetime will become self-similar only after the passage of a sufficent 
amount of time during which the original differences will be washed away and the collapse 
approaches a self-similar evolution, such that the physics on large scales is related to that on 
small scales in a universal manner. A schematic illustration of this process is given in Fig. 
8. Region A represents a region where the influence of the terms V(4>) and A in Einstein 
field equations is still very strong, and the corresponding spacetime is significantly different 
from that of the critical solution. In particular, no trapped surfaces have been formed so 
far. As r grows, the differences are smoothed out, and when r grows to a certain value, 
say, t = Ti, the differences become very small and can be considered as perturbations on 
top of the the self-similar critical solution. The region where such perturbations are valid 
is denoted as region B in Fig. 8. Clearly, the time T\ is dependent on the initial-data, and 
by fine tuning one can find out which initial data leads to critical collapse and which does 
not. When the strength of the initial data is larger than that required for critical collapse, 
the collapse will be strong enough to form black hole, while when weaker, the scalar field 
will reflect through the origin before a black hole can be formed and subsequently disperse 
entirely to infinity. In region B, the solutions to the whole problem can be expanded in the 
form PI, 
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(150) 



n=0 



where A = (a, s, tp), and A* denotes the critical solution. Substituting the above expres- 
sion into the Einstein field equations (|TT|-|2"T|) , and then separating terms by powers of e~ r , 



we find that, to zeroth order, the resultant equations are exactly those given by Eqs. (27-31) 



for the self-similar critical solution A^\ which are independent of other higher terms. On 
the other hand, the equations for A* s with n > 1 are always coupled to its lower terms, 



A* , ■ ■ ■ ,A ( *'. This allows us to determine A* l) recursively starting with A^' . Thus, m 
principle once A^ is known, we can find all the higher terms A*. In particular, one can 
show that A* satisfies the same equations as those for A\{z) given by Eq. ([77H81]) with 
k = —1, where A\(z) = (s\, a%, <£>i). Following the analysis given there, we find that 
A^ is well-behaved in region B, regular at r = 0, and analytical across the hypersurface 
v = 0. Since the higher order terms A^ satisfy similar types of differential equations as 
those for A^\ and also these higher terms in ( |150| ) all correspond to stable modes (k < 0) 
of perturbations given by Eq.([T|), it is quite reasonable to expect that all the higher order 
terms in Eq. (|150|) are well behaved in Region B, and in the limit r — > oo Eq.(150) converges 
to A^ . Note that A^ is well-behaved in the whole region B, except on the hypersurface 
v = 0. However, this divergence is weak (logarithmically) and is due entirely to the choice 
of coordinates. As shown in section III, in terms of u and v the metric cofficients are well 
behaved there. In review of all the above, we can see that Eq.( |150| ) indeed represents solu- 
tions to the full problem, that is, solutions of Einstein-scalar equations with potential V^(^) 
and the cosmological constant A. 



|(0) 



i(n) 



|(0) 
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FIG. 8. Global structure of near-critical gravitational collapse. In region A the effects of the 
terms V(4>) and A are strong and the spacetime will, in general be very different from that of the 
critical solution. In particular, an apparent horizon will not have formed. With increasing time, 
however, the effects of the potential will become more and more weak such that, in region B and 
for those initial configurations with initial data close to criticality, those effects can be considered 
as perturbations of the critical solution. In this region, supercritical configurations will form, in 
the end, black holes, while subcritical configurations will reflect through the origin and disperse 
to infinity without formating black holes. The line AH represents the apparent horizon formed in 
the supercritical collapse. 



Now let us consider the linear perturbations of the solutions ( |150| ), which can be written 

as 

5A(r,z) = e kT 5 k A(T,z), (151) 

where SkA(r, z) denotes the general perturbations of A(t, z) for the mode k, similar to that 
given by Eq.(|75|), but now depending on both of r and z. However, we can expand it in a 
similar form as that given by Eq.( |15(J| ), 

CO 

5 k A{r,z) = Y,e- nT 5 k A^\z). (152) 

n=0 

As A(t, z) is completely determined by A^\ 8kA(r, z) is totally specified by its leading 
term SkA^(z). Indeed, substituting Eqs.( |150| - |152| ) into the Einstein field equations (|Hj- 



21"D , we find that S k A^(z) satisfies the same differential equations flTTHBlD for Ai(z), where 
the background depedence is only on terms of A^\ Consequently, the spectrum for the 
perturbations of S k A^(z) is the same as that of Ai(z). Once 5 k A^(z) is known we can 
find the higher order terms 5 k A^ n \z) recursively. In this way, we can see that the spectras 
of perturbations in the cases with or without V((p) and A will be exactly the same. In 
particular, the critical solution for the massless scalar field is also a critical solution for 
a scalar field with potential in a non-vanishing cosmological constant background (with the 
possibilities of being either negative or positive), provided Eg. holds. 



VI. CONCLUSIONS 

In this paper, we have considered the Einstein field equations for a scalar field with 
potential V(<p) in a background where the cosmological constant A is different from zero. 
We have argued that V(4>) and A will not affect the critical behavior provided the condition 
(f26|) is true. 

All the continuously self-similar solutions to Einstein-massless-scalar field equations were 
found together with their analytic extensions. It was shown that some of these solutions 
represent gravitational collapse of the scalar field. In all these collapsing models, a black hole 
is always formed. It should be noted that recently Ida showed that no black holes can be 
formed in (2+1) dimensional gravity, unless the dominant energy condition is violated |L5| . 



This seems to conflict with the results presented in this paper. One possible explanation 
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may be that the conditions assumed in [[HJ do not hold here, specifically the one that the 
event horizon is a compact surface. 

In section IV, we also studied the linear perturbations of all the self-similar solutions dis- 
cussed in section III, and found all the perturbations in closed form. On imposing boundary 
conditions of (1) regularity at the origin and (2) analyticity and an inflow only condition at 
the event horizon, we were able to identify the unique, exact solution which has a single un- 
stable mode. This solution is that given by Eqs.flo7|) and flBBQ with n = 4, and which was first 
found by Garfinkle JF|, and which he showed matches well with the critical solution found 
numerically by Pretorius and Choptuik While this result is certainly very encouraging, 
we must note that in our case the exponent 7 that can be obtained from the expression 
7 = 1/k — 4, is significantly different from the one determined numerically by Pretorius and 
Choptuik, 7 ~ 1.2 ±0.05. Nor is it similar to the critical exponent of 7 ~ 0.81 which Husain 
and Olivier || find using a different numerical code written in double null coordinates. The 
recent perturbative results of Garfinkle and Gundlach ||, while finding the n = 2 solution 
to be the critical solution, do get a better critical exponent, namely 7 = 4/3 ~ 1.33. All of 
these differences seem quite intriguing and together with the problem as a whole, are worthy 
of further investigation. 
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APPENDIX: EXACT SOLUTIONS FOR MASSLESS SCALAR FIELD 

As shown in section II, when r — ► 00, the terms associated with the potential V(4>) always 
go like e~ 2T V((p). Thus, for the cases where V(4>) is finite, the potential terms are negligible 
in the strong field regime, where critical behavior becomes visible. Setting V(4>) = in 
Eqs.(|T3Hl5|), we find that the corresponding Einstein-scalar field equations become, 

r >uu - 2a >u r >u = -8nGr(p 2 u , (A.l) 

r> - 2o>r> = -8rrGr(p 2 v , (A.2) 

r )UV + 2ra tUV = -87rGr<p )U <f> >v , (A. 3) 

r, uv = 0, (A.4) 

+ - {r,u<P,v + r >v (f), u ) = 0. (A.5) 



,uv 

r 



From Eq.( |A.4j ) we find that the metric coefficient r has the general solution, 

r = a{u) + (3{v), (A.6) 

where a(u) and f3(v) are arbitrary functions of their indicated arguments. Clearly, in general, 
the solutions can be divided into three different cases, 
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(i) a'(u) = 0, (3'(v) ^ 0; (ii) a'(u)/3'(v) ^ 0; (Hi) a'(u) ^ 0, f3'(v) = 0. (A.7) 



In the cases (i) and (in), by exchanging the two variables u and v, we can get one from the 
other. Therefore, without loss of generality, in the following we shall consider only Cases (i) 
and (ii). 



A.l Exact Solutions for ot!(u) = 0, (3'{v) ^ 



In this case we have a(u) = const. By redefining the function /3(v), we can always set 
this constant to zero. Then, from Eqs. flA.ll ) and (|A.3|) we find that 

= 0(ti), a = a(u) + b(v), (A. 8) 

where a(u) and b(v) are other two arbitrary functions. However, by introducing new coor- 
dinates u and v via the relations, 



du = e 2a(u W dv = e 2b(v) dv, 



(A.9) 



we can always set a(u) = = b(v), a choice that will be assumed in the following discussion. 
Then, substituting Eq.( [A.8|) into Eq.( [A.2|) , we find that 



<f>(v) 



V8ttG 



(3(v) 



1/2 



dv + 4> l Q . 



(A.10) 



Thus, for any given function /3(f), Eq. flA.lOD completely determines the massless scalar field, 
provided that 



(3(v) > 0, (3"(v) < 0. 



(A.ll) 



On the other hand, in order not to have the spacetime be closed at the beginning of collapse, 
we also require that 



(3'(v) > 0, (v > v Q ), 



(A.12) 



for some value of Vq, where Vq is a constant. Combining Eqs.( A.ll ) and ( |A.12 ), it can be 
shown that in the present case the collapse never forms black holes. This is because the 
formation of black hole is indicated by the formation of an apparent horizon, on which we 
have 



P'{vah) = 0. 



(A. 13) 



However, Eqs. flA.ll^ ) and ( |A.13| ) with the condition that f3(v) > imply that f3"(v) > 
right outside the apparent horizon. The latter is possible only when is imaginary, which 
is forbidden by the weak energy condition [[16]. 

It should be noted that, although the collapse cannot form black holes, it can form 
singularities that are null and strong. To show this, let us first notice that in the present case 
all the scalars built from the Riemann curvature tensor are zero, therefore, scalar curvature 
singularities are always absent [17]. However, we do have non-scalar curvature singularities 
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at r = (or equivalently (3{y) = 0), in the sense that the tidal forces experienced by an 
observer become infinitely large as the hypersurface r = approaches. To this end, let us 
consider the radial timelike geodesies, which in the present case are simply given by 



u 



,—2do 



= 0, 



(A.14) 



where do is an integration constant, and an overdot denotes the ordinary derivative with 
respect to the proper time, A, of the timelike geodesies. Defining ej™ = dx a /d\, we find that 
the three unit vectors, given by 



'(2) 



-2d ra 



1 



-2d ra _ i p 2d ra 

-On, 



(A.15) 



form a freely falling frame, 



e 1c) e \d)9ab 



//,,/. <: ::l/ ': iU 0, (A.16) 

diag. { — 1, 1, !}• Projecting the Ricci tensor onto the above frame, we find 



where r] ab 
that 



R 



(a)(6) 



Rr 



(a) e (6) 



4d 



(A.17) 



To study the asymptotic behavior of the solutions near the hypersurface r = (3{v) = 0, 
we expand the function (3{v) near this surface and assume that the leading term of the 
expansion is proportional to t> 7 , i.e., 



(A.18) 



where the constants (3q and 7 have to satisfy the conditions (3q > and < 7 < 1, in order 
to have Eqs.( [TOl| ) and ([A~T^ ) hold. Substituting Eq.( |AT8D into Eq. flATTTD , we find that 



R («)W = 1757(1 - 7)e 4d0 [6°J° b - (5°Jl + + 6161 



Av 
7(1-7) 
A 2 



SnSh - ( 8*8h + Oh ) + <^M, 



'a u b 



cOrl 



IrO 



J a u b 



(A.19) 



where the proper time A along the timelike geodesies was chosen such that r = corresponds 
to A = 0. The above equation shows clearly that R( a )(b) becomes unbounded as r — > 0. Thus, 
the hypersurface r = represents a non-scalar spacetime singularity. This singularity is 
strong fll8| , in the sense that the distortion, which is proportional to two integrals of R( a )(b) 
with respect to the proper time A, becomes unbounded 



dX I i?( a )( b )C?A 



7(7-l)m(A) 



00, 



(A.20) 
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as A — > 0. Note that this is different from what was claimed in [|T^ . 

Before turning to the next case, let us note that the self-similar solutions given by Eq.i 
correspond to the choice, 



(3{v) 



v 



V(x+i) 



(A.21) 



A. 2 Exact Solutions for a'(u)(3'{v) / 



In this case solving Eqs.( [A.i] - |A.5| ) in general is very difficult. Instead, we can prove the 
following theorems. 

Theorem 1: For any given solution, M(u,v), U(u,v), and V(u, v) of the Einstein 
vacuum field equations, Rffl = 0, for the four-dimensional spacetimes described by the 
metric, 



aS (4) 



the solutions, 



a = — 

r = e~ 



1 r 

2 

■u 



-2e- M dudv + e~ u (e v dx 2 + e~ v dy : 



(2 - 3 7 2 )£/ - (1 - 2 7 2 ) In \U U U V \ + 2 7 2 M 



(A.22) 



o"o, 



7 



V. 



(A.23) 



satisfy the Einstein-scalar field equations ( |A.1| - |A.5| ), where 7 and Uq are two arbitrary con- 
stants. 

The proof of the above theorem is straightforward. In fact, substituting Eq.( [A.23"|) into 
Eqs.( A.l - |A.5|) , we find that these equations, in terms of M, U and V, take the form 



2U UU - U% + 2U U M, U = V 2 U , 
2U VV - U 2 V + 2U V M >V = V 2 , 
2M tUV + 2U UV - U, U U, V = V U V V , 
2V UV - U U V V - U V V U = 0, 



(A.24) 
(A.25) 
(A.26) 
(A.27) 
(A.28) 



which are exactly the Einstein vacuum field equations for the four-dimensional spacetimes 
described by the metric QA.22 ) [19 . Introducing a new timelike coordinate t via the relation 



t = a{u) — P(v), 



(A.29) 



we find that t and r are linearly independent, because now we have the condition a'(u)f3'(v ) 7^ 
0. Then, in terms of t and r, the metric ( |A.22| ) can be written in the form 



d? 2 

aS (4) 



3 n(t,r) 



- dr 2 ) + r [e v ^dx 2 + e~ v ^ r W 



(A.30) 



where 
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Q = —M — -U — In \2a'(u)P'(v)\ . 



Hence, in terms of Q the Einstein vacuum field equations ( |A.24| - [A~.28|) take the form 

= \ {v 2 r + vf) , 



n 



n 



r 
2 

rV,Vu 



a, 



2V rr + -V r - V tt = 0. 

r ' 



(A.31) 

(A.32) 
(A.33) 
(A.34) 

(A.35) 



On the other hand, substituting Eq.( |A.3l| ) into Eq.( |A.23| ), we find that 

a = 1 2 Q + hn\a'(u)(3'(v)\+a , 
t = a(u) — f3(v), r = a(u) + (3(v), 



7 



V, 



(A.36) 



By properly imposing physical conditions J20|; the metric given by Eq. (|A.22| ) can be inter- 
preted as representing cylindrically symmetric spacetimes with one of the two coordinates 
x and y being the angular coordinate and the other being the axial coordinate. But, if we 
take the values of x, y in the range — oo < x, y < oo, the metric is usually considered 
as representing cosmological models ||21|| . However, since here we are mainly interested in 
obtaining exact solutions of the (2 + l)-dimensional metric (|5]), we shall not be concerned 
with the physical interpretation of the metric ( |A.22j ). 

It is interesting to note that if we make the replacement 



in Eqs.( [A~32HA^ ) we find that 

n z = 



t = iz, 



vi-v: 



rV r V z , 



2V„ + -V r + V z 



0. 



(A.37) 

(A.38) 
(A.39) 
(A.40) 

(A.41) 



which are exactly the Einstein vacuum equations for four- dimensional spacetimes with ax- 
isymmetry |2"T| 



e ^-v (d r i + dz ^ + r 2 e- y d^ 2 - e v dt 2 . 



(A.42) 



Theorem 2 [2^j: For any given solutions M(u,v), V(u,v), and U(u,v) of the Einstein 
vacuum field equations ( |A.24bfA~28|) , those with 
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F(u, v) = (l- 5 2 ) [-U-ln \U U U V \ )+5 2 M + F , 



V. 



(A.43) 



where 5 and Fo are two arbitrary constants, satisfy the Einstein-scalar field equations, R$ = 
8^(7$^$^, of the four- dimensional spacetimes with plannar symmetry, 

(A.44) 



ds\ 4) = -2e~ F ^ v) dudv + e~ u(u ' v) (dx 2 + dy 1 
The combination of Eqs.( [A.2"3| ) and ( |A.43| ) yields, 

^ {(3 7 2 - 25 2 )U + (5 2 - 2 7 2 ) In \U U U V \} - ^F 
V2~l 



r[u,v) 



25 2 



o- 



[U,V 



(A.45) 



Since the vacuum spacetimes of Eqs.( |A.2"2"D , ( |A.42| ) and the massless scalar spacetimes 
of Eq.( [A.44|) have been extensively studied, and many analytical solutions are known 
|22] , |2l|j23| , [24]] , we can use the above two theorems to construct solutions of the Einstein- 
scalar field equations in (2 + l)-dimensional spacetimes described by the metric (|5|). In 
particular, Eq. (|A.27|) has the general solutions, 



2 N N 

V = ct + d ln(r) + - ^ g k In (a k ) + ^ /i fc Arcos 

^ k=l k=l 



2t k a k 



1/2 



+ 



+ 



r(l + a k ) 

[A(uj) sin(a>£) + B{oj) cos(a;t)] Jo {uir)duj 
\C{tjj) sin(u;£) + D(uj) cos(a;t)] N (ur)du, 



(A.46) 



where c, d, g k , h k , A(u), B(u), C(u), D(u) and uj are arbitrary real constants, and N is 
an integer. The functions Jo {ujr) and N {ojr) denote the Bessel and Neumann functions of 
zero order, respectively, while a k can be chosen as o~ k or o~ k , which are defined by 



= L k ± (L\ - l) , 



+ 



1 - 



(A.47) 



where 



t k = t° k -t, (k = 1, 2, N), (A.48) 
with t k and w k being real constants. From the above expressions it can be shown that 

< a k < 1, 1 < a£ < oo. (A.49) 
The third and fourth terms in Eq.( |A.46| ) are usually called generalized soliton solutions. 
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The corresponding general solutions for Q are not known, but, in some particular cases 
one can find it from Eqs.( A.38| ) and ( |A.39| ) by quadratures. For example, when d and g k s 
are the only constants that are different from zero, it is given by, 



for 



(d 2 ~g 2 )/2 



N 

n 



a 



g k (2g k +d-g)/2 



N 



k=i [(1 - a k ) 2 H k ] 9k/4: ijJi^j) 



n 



gtgj/2 



(A.50) 



where f is a constant, and 



A' 



9 = J2dk, H k = (1 - a k f + — — 



k=l 



Ok) 



Hj 



(<Ji + aAr 



ItitjCriCTj 



-i 2 



[l + ai )(l + a 



3)] 



(A.51) 



(l + «7i)(l + 0- i;j 

On the other hand, if d and h k s are the only non-zero constants in Eq. ([A.46|) , the 
corresponding Q is given by 



for 



N 



(d 2 +2h)/2 



N 



n # 

M=l(i>j) V 



\ hihj/2 



N 

x exp < d ^ /i^Arcos 
fe=i 



2t fc (7 



1/2 



r(l + <7 fc ) 



2^ fe (fe fc -l)/4 



k=i n k 



h k (h k -2)/A 



r 2 (l - a k )] n - 



hl/2- 



where H k is defined by Eq. ( |A.5ll ), and 



N 



h = Y,h k , A k = (t 2 k -w\- r 2 ) 2 + Aw 2 t 



2+2 

ki 



k=l 



Af^{l-a 2 ){l-a 2 ){r 2 {a l + a 3 ) 



tjtj 



± 



WiWj 



_(l + <7i)(l + ^) (l-ai)(l-<Tj)\j' 
Moreover, the separation of variables of Eq. (|A.41|) yields the solutions [p4 



jr I \ -Pn(cOS 

V A xiai{r, z) = 2^ a n - 



n=0 



+ 1 



where (Xn £1X6 constants, and 



_ 2 i 2 

p = r + z , 



cos# = -. 

p 



Then, the integration of Eqs.( |A.35D and ( |A.39| ) gives 
^Axiai(r,z) 



E a n a lBln , p p p p 



n «+n+2 
l,n=0 " 

where P n represents the Legendre polynomial of order n, and 

(Z + l)(n + l) 



An = 



2(Z + n + 2) 



(A.52) 



(A.53) 



(A.54) 



(A.55) 



(A.56) 



(A.57) 
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